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Abstract 

A theory of spontaneous parametric down-conversion, which gives rise to 
a quantum state that is entangled in muhiple parameters, such as three- 
dimensional wavevector and polarization, allows us to understand the unusual 
characteristics of fourth-order quantum interference in many experiments, in- 
cluding ultrafast type-II parametric down-conversion, the specific example 
illustrated in this paper. The comprehensive approach provided here permits 
the engineering of quantum states suitable for quantum information schemes 

and new quantum technologies. 
42.50.Dv, 42.65.Re, 42.65.Ky, 03.67a 
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Entanglement |l| is, undoubtedly, one of the most fascinating features of quantum me- 
chanics. Spontaneous parametric down-conversion (SPDC) a nonlinear optical phe- 
nomenon, has been one of the most widely used sources of entangled quantum states. In 
this process, pairs of photons are generated in a state that can be entangled in frequency, 
momentum, and polarization when a laser beam illuminates a nonlinear optical crystal. The 
experimental arrangement for producing entangled photon pairs is simple both in conception 
and in execution. 

Ironically, a significant number of experimental efforts designed to verify the nonsepara- 
bility of entangled states, the hallmark of entanglement, are carried out in the context of 
models that fail to access the overall relevant Hilbert space, but rather are restricted to only 
a single kind of entanglement, such as entanglement in energy 0], momentum or polar- 
ization [Q. Inconsistencies in the analysis of down-conversion quantum-interference experi- 
ments can emerge under such circumstances, as highlighted by the failure of the conventional 
theory of ultrafast parametric down-conversion to characterize quantum-interference ex- 
periments 101 . 

In this paper we present a quantum-mechanical analysis of entangled-photon state 
generation via type-II SPDC, considering simultaneous entanglement in three-dimensional 
wavevector and polarization at the generation, propagation, and detection stages. As one 
specific example of the applicability of this approach, we use it to describe both new and 
previously obtained 0] results of SPDC experiments with a femtosecond pump. Our analysis 
confirms that the inconsistencies between existing theoretical models and the observed data 
in femtosecond down-conversion experiments can indeed be attributed to a failure of con- 
sidering the full Hilbert space spanned by the simultaneously entangled quantum variables. 
Femtosecond SPDC models have heretofore ignored transverse wavevector components and 
have thereby not accounted for the previously demonstrated angular spread ||^ of the down- 
converted light. The approach presented here is suitable for Type-I, as well as Type-II, 
spontaneous parametric down- conversion in the paraxial approximation, which is valid for 
the great preponderance of experimental SPDC efforts to date. 
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Our study leads to a deeper physical understanding of hyperentangled photon states 
and, concomitantly, provides a route for engineering these states for specific applications, 
including quantum information processing. 

Hyperentangled- State Generation. — With this motivation we present a multidimensional 
analysis of the entangled-photon state generated via SPDC To admit a broad range of 
possible experimental schemes we consider, in turn, three general and fundamentally distinct 
stages in any experimental apparatus: the generation, propagation, and detection of the 
quantum state 

We begin with generation. By virtue of the weak nonlinear interaction, we consider the 

state generated within the confines of first-order time-dependent perturbation theory: 

t 

1^^'^) ~ j dt' H,^,{t') |0). (1) 

to 

Here Hint{t') is the interaction Hamiltonian, {to,t) is the duration of the interaction, and |0) 
is the initial vacuum state. The interaction Hamiltonian governing this phenomenon is 

^int(t')~X^'^ / rfrEW(t', r)Ei-\t', v) E^-\t' , r) + H.c. , (2) 

V 

where x''^^ is the second-order susceptibility and V is the volume of the nonlinear medium 
in which the interaction takes place. The symbol Ej'^\t', r) represents the positive- 
(negative-) frequency portion of the jth electric-field operator, with the subscript j repre- 
senting the pump {p), ordinary (o), and extraordinary (e) waves at time t' and position 
r, and H.c. stands for Hermitian conjugate. In the paraxial approximation, the polariza- 
tion of each photon (o, e) may be assumed to be independent of frequency and wavevector 
(the dependence at large angles is considered in Ref. Because of the high intensity 

of the pump field we take the coherent-state laser beam to be classical, with an arbitrary 
spatiotemporal profile given by 

Ep{v, t)= I dk, ^p(kp)e'''--'-e-''^''('^'')* , (3) 

where £'p(kp) is the complex-amplitude profile of the field as a function of the wavevector 



In contrast with previous models we consider the wavevector to be three-dimensional, 
with a transverse wavevector qp and frequency cUp, so that Eq. (|^) takes the form 

Ep{v, t)= j dqp dup Ep{qp, u;p)e'''-'e'''--''e-'''-' , (4) 

where x spans the transverse plane perpendicular to the propagation direction z. In a similar 
way the signal and idler fields can be expressed in terms of the quantum-mechanical creation 
operators a^(q, (j;) for the (q, cj) modes as 

where the subscript j = o, e. The longitudinal component of k, denoted k, can be written 
in terms of the (q, u) pair as |^ 



\ 



, (6) 



where 6 is the angle between k and the optical axis of the nonlinear crystal, ne{uj,9) is 
the extraordinary index of refraction in the nonlinear medium, and c is the speed of light 
in vacuum. Note that the extraordinary refractive index, ne{u,6), in Eq. (j^) should be 
replaced by the ordinary refractive index, no{uj), when calculating k for ordinary waves. 

Substituting Eqs. (||) and (|) into Eqs. (|l]) and (||) yields the wavefunction at the output 
of the nonlinear crystal: 

|^(^)) ~ J dqodqeduJodiOe ^{cio,cie]^o,^e)dl{qo,u;o)al{qe,uje)\0) , (7) 

with 

<I>(qo,qe;a;o,u;e) = Ep{qo + qe](^o + ^e) Lsmc{^)e~'^ . (8) 

Here A = where Hj {j = p,o,e) is related to the indices [qjjUj) via relations 

similar to Eq. (||). The nonseparability of the function $(qo, qe; i^o, i^e) in Eqs. (|^) and (||), 
recalling (|^), is the hallmark of simultaneous multi-parameter entanglement. 

Hyperentangled- State Propagation. — Propagation between the planes of generation and 
detection is characterized by the classical transfer function of the optical system. The 



biphoton probability amplitude at the space-time coordinates (xyi,tyi) and (x^jts), where 
detection will take place, is defined by |I0| , 



The explicit forms of the quantum fields present at the detection locations are represented 

by 

E^/\xA,tA) = J dqduj e""^*-^ [7^Ae(xA, q; cu)ae(q, + 7^ao(xa, q; ^)ao(q, tc')] , 
E^j^\-i^B,tB) = j df\duj e""^*^ [7^Be(xiJ,q; a;)ae(q,u;) +7^Boi?^B,^:,^)do{,^.,^^)] , (10) 

where the transfer function Tiij {i = A,B and j = e, o) describes the propagation of a (q, lo) 
mode from the nonlinear-crystal output plane to the detection plane. Substituting Eqs. (|^ 
and dTDI ) into Eq. (|^) yields a general form for the biphoton probability amplitude: 

A(xa, xb; tA, Ib) = J dqodqe duodue $(qo, qe; ^o, tUe) 

X HAei^A, qe; UJe)nBo{^B, ^o] ^o) e-'("=*^+"°*^) 

+ Hao{^A, qo; a;„)7^Be(xB, qe; UJe) e-'(-°*A+..eiB)] _ (11) 

This function can be separated into polarization-dependent and -independent terms, as 
necessary, for any particular configuration. By choosing explicit forms of the functions TiAei 
'Hao, 'Hbc, and Hbo, the overall biphoton probability amplitude can be sculpted as desired. 

Hyperentangled- State Detection. — The formulation of the detection process depends on 
the scheme to be used. Slow detectors, for example, impart temporal integration while finite 
area detectors impart spatial integration. Quantum-interference experiments typically make 
use of just such detectors. Under these conditions, the coincidence count rate R is readily 
expressed in terms of the biphoton probability amplitude: 

R = J dxAd^BdtAdtB |A(x^,XB;tA,^i?)r • (12) 

Example: Quantum Interference in Ultrafast SPDC. — We now consider a particular ex- 
ample that demonstrates the validity of our analysis: an ultrafast polarization quantum- 
interference experiment of the form illustrated in Fig. 1(a). Details of the experimental 



arrangement and protocol can be found in an earlier work [0; in the analysis offered there 
we made use of a phenomenological model that considered a collection of contributions from 
different regions in the nonlinear crystal that, in the absence of a full quantum- mechanical 
model, were conjectured to be independent and distinguishable. With the help of the general 
spatiotemporal quantum-mechanical approach developed here, we are now in a position to 
provide a complete analysis of those data along with new data in which filtering was used, 
presented in Figs. 2 and 3, respectively. 

For the polarization-interferometer arrangement illustrated in Fig. 1(a), in the pres- 
ence of a polarization-dependent relative temporal delay r, Eq. (|l^) can be conveniently 
separated into polarization-dependent and -independent terms via the relation 

7^,,(x„q;u;) = (e, ■ e,) e"'"^^- i/(x„ q; a;) , (13) 

where i = A,B and j = e, o. The symbol 6ej is the Kronecker delta so that See = 1 
and Seo = 0. The unit vector ej describes the orientation of the polarization analyzers in 
the experimental apparatus [see Fig. 1(a)], while ej is the unit vector that describes the 
polarization of the down-converted photons; the function H{:Ki, q; u) is the transfer function 
of the polarization-independent elements of the system such as free space, filters, apertures, 
and lenses, as illustrated in Fig. 1(b). The paraxial approximation satisfactorily describes 
our experiments so that the explicit form of H in Fig. 1(b) becomes 



cf 



H^) , (14) 



where di, (i2, and / (focal length of the lens) are indicated, P is the aperture function p{x) 
in the Fourier domain, and J^{uj) is the spectral filter function. 

Using Eqs. (p!3|) and (|T^) in Eq. (pH]), the biphoton probability amplitude for the 



arrangement shown in Fig. 1(a) therefore becomes 

v4(xa, Xb; Ia, ts) = J dqodqe dujodue $(qo, qe] oJo, cJe)e~"^"^ 

X "(e^ ■ ee){eB ■ e„) H{^a, qe; oJe)H{^B, qo; u;o)e~'("^*^+"°*^) 

+ (e^ ■ eo)ieB ■ e^) //(x^, qo; uJo)HixB, qe; a;e)e-'("°*-^+"^*^)l . (15) 



Using this form for the biphoton probabihty amphtude in Eq. (JT^) yields the coincidence- 
count rate as a function of the polarization-dependent temporal delay r. 

Discussion. — Figure 2 displays the observed normalized coincidence rates (fourth-order 
quantum-interference patterns) for 0.5-, 1.5-, and 3.0-mm BBO crystals (symbols), in the 
absence of spectral filtering, along with the expected theoretical curves (solid), as a function 
of relative optical-path delay r. We have treated the pump as a finite-bandwidth pulsed 
plane-wave, an assumption that is valid in our experimental setup. The asymmetry of the 
observed interference pattern clearly increases with crystal thickness. 

Figure 3 provides a set of data collected in a similar fashion, but this time observed in 
the presence of a narrowband (9-nm) spectral filter J^{uj), as illustrated in Fig. 1(b). The 
most dramatic effect of including the filter is the symmetrization of the quantum-interference 
patterns. Since q and u are intrinsically linked by Eq. (H), the imposition of spectral filtering 
restricts the allowable transverse wavevector spread. Spectral and spatial filtering therefore 
have similar effects for non-cross-spectrally pure light, such as that generated in SPDC [Q. 

The increasing asymmetry and loss of visibility in Figs. 2 and 3 are observed with increas- 
ing crystal thickness, as the extent of the (q, u) modes overlap less in space at the detection 
plane. This decreased overlap leads to increased distinguishability. This distinguishability 
is similar in nature to the spectral distinguishability in the one- dimensional model discussed 
in Ref. 0. The physical origin of this behavior resides in the angle- dependence (hence q- 
dependence) of the extraordinary refractive index for the down-converted photons [Eq. (|^)]. 
Although the phase-matching condition between the pump and the down-converted photon 
pairs encompasses a large range of (q, u) modes at the source, the combination of free- 
space propagation and the small acceptance angle of the optical system leads to diffraction 
of the SPDC beams, which, in turn, results in increased overlap and therefore a decrease 
in distinguishability. Indeed, when the aperture size becomes sufficiently small, the ob- 
served quantum-interference patterns ultimately revert to those calculated using the one- 
dimensional model that has traditionally been employed. 
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FIGURES 

FIG. 1. (a) Schematic of the experimental setup for observation of quantum interference using 
femtosecond SPDC. (b) Detail of the path from the crystal output plane to the detector input 
plane. 

FIG. 2. Experimental (symbols) and theoretical (solid curves) results for the normalized coinci- 
dence rate for BBO crystals of three different lengths (hexagons: 0.5 mm; triangles: 1.5 mm; circles: 
3.0 mm) as a function of the relative optical-path delay r. As the crystal length increases the fringe 
visibility diminishes substantially and a dramatic asymmetry emerges. No free parameters are used 
to fit the data. 

FIG. 3. Plots similar to those in Fig. 2 in the presence of an interference filter of 9-nm band- 
width. The patterns are symmetrized. 
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